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1.  Introduction 


It  has  been  shown  that  the  efficacy  of  radiotherapy,  chemotherapy  and 
photodynamic  therapy  are  influenced  by  tumor  oxygenation,  and  hypoxic  tumors  are 
refractory.  Thus,  a  reliable  non-invasive  method  of  measuring  tumor  oxygen  tension 
would  be  valuable  in  treatment  planning.  Many  techniques  for  monitoring  pO^  have 
been  developed.  Most,  however,  are  invasive  {e.g.,  electrodes,  biopsy  and  cryospec- 
trophotometry)  or  limited  to  surface  tissues  (fluorescence  and  phosphorescence). 
MRI  of  hexafluorobenzene  (HFB),  based  on  Pulse  Burst  Saturation  Recovery  (PBSR) 
-  Echo  Planar  Imaging  (EPI),  provides  an  alternative  non-invasive  technique  to 
monitor  pO^  or  determine  the  extent  of  hypoxia.  Implementation  of  this  technique, 
however,  requires  a  high  quality  double-tuned  (’H  and  '®F)  whole  body  birdcage 
resonator  or  a  continuously  swithable  whole  body  resonator  in  order  to  niiiimiize 
preparation  time  and  prevent  accidental  death  of  rat.  The  desired  imaging  and 
spectroscopy  resonator  for  our  purposes  would  be  the  one  that  provides  the  highest 
possible  signal-to-noise  ratio  (SNR),  a  good  filling  factor,  and  a  highly  homogeneous 
Bi  field.  To  address  these  issues,  I  designed,  built,  and  tested  two  RF  resonators,  one 
being  a  double-tuned  whole  body  birdcage  resonator  and  another  one  a  swithable 
whole  body  slotted  tube  resonator. 

For  the  birdcage  resonator,  laboratory  bench  testing  indicated  that  the  unloaded 
Qs  for  both  proton  and  fluorine  resonant  modes  are  similar  and  relatively  low.  The 
loaded  Q  drops  by  about  36%  for  fluorine  and  39%  for  proton.  Phantom  imaging 
found  that  the  Bi  fields  are  somewhat  heterogeneous  across  usable  volume  of  the 
resonator  for  both  resonant  modes.  For  the  slotted  tube  resonator,  laboratory  bench 
testing  shown  that  the  resonator  has  the  capability  of  being  continuously  tuned  in  the 
frequency  range  of  150  MHz  ~  220  MHz.  The  resonator  has  high  Q  values  and  short 
90°  pulse  widths  for  both  proton  and  fluorine,  which  results  in  improved  sensitivity 
and  makes  in  vivo  tumor  pO^  studies  possible.  In  addition,  the  Bt  field  homogeneity 
was  found  to  be  excellent  for  proton  resonance  and  reasonably  good  for  fluorine 
resonance  across  usable  volume  of  the  resonator. 

HFB  has  been  found  to  clear  relatively  rapidly.  Clearance  rate  may  provide  an 
indication  of  tumor  perfusion.  To  compute  the  clearance  rate  on  a  voxel  by  voxel  basis 
and  thus  produce  a  clearance  map,  I  wrote  a  Unix  based  NMR-shell  script  program.  In 
addition,  I  wrote  another  NMR-shell  script  program  to  display  the  T1  relaxation  curve 
and  assess  the  goodness  of  the  curve  fit  of  individual  voxels  of  PBSR-EPI  images.  The 
graphical  representation  of  the  data,  clearance  model,  relaxation  model,  and  goodness 
of  fit  provided  a  quick  and  convenient  way  of  assessing  the  quality  of  the  data.  Both 
programs  have  been  tested. 


2.  RF  Resonator  Design  and  Testing 


2.1  Functions  and  Characteristics  of  RF  Resonators 

The  fundamental  functions  of  an  NMR  radiofirequency  (RF)  resonator  or  coil 
are  two  fold:  (1)  to  excite  the  nuclear  spins  of  interest  by  transmitting  a  magnetic  field 
Bi  and  then  after  a  suitable  delay,  (2)  to  receive  the  resulting  signal  from  the 
precessing  nuclear  spin  magnetization,  called  free  induction  decay  (FID).  During 
excitation  mode,  the  RF  resonator  serves  as  a  transducer  that  converts  RF  power  into 
the  transverse  rotating  magnetic  field  Bi  in  the  imaging  volume  of  a  sample. 
Therefore,  high-transmission  efficiency  means  that  the  resonator  delivers  a  maximum 
Bi  to  the  sample  volume  with  a  minimum  RF  power.  During  reception  mode,  the  RF 
resonator,  along  with  its  associated  preamplifier,  acts  as  a  transducer  that  converts  the 
precessing  nuclear  spin  magnetization  into  electrical  signals  suitable  for  further 
processing  by  a  computer.  High-detection  efficiency  for  the  reception  mode  then 
corresponds  to  minimal  degradation  of  the  inherent  signal-to-noise  (SNR)  of  the 
sample  volume.  A  well-designed  and  constructed  RF  resonator  should  be  highly 
efficient  as  a  transmitter  as  well  as  a  receiver. 

The  NMR  signal  (FID)  is  a  time-dependent  AC  voltage  induced  in  the  RF 
resonator  by  the  nuclear  spin  magnetization  vector  M  precessing  freely  in  the 
transverse  (x’y’)  plane.  The  M  has  been  nutated  into  the  jc’y’  plane  by  a  torque  exerted 
on  it  by  the  rotating  (pulsed)  magnetic  field  jB,  of  appropriate  frequency  (Larmor 
frequency),  time  duration,  and  amplitude.  The  RF  resonator  transmits  an  RF  power  of 
the  order  of  hundreds  of  volts  into  the  sample  volume,  whereas  it  receives  a  voltage 
from  the  magnetization  M  of  the  order  of  one  tenth  of  a  microvolt  (0. 1  pV). 

The  relative  orientation  of  the  static  magnetic  field  Bo,  the  rotating  RF  magnetic 
field  Bi  and  the  precessing  magnetization  vector  M  is  crucial  in  NMR  experiments. 
The  Bi  filed  must  be  perpendicular  to  Bo,  and  the  RF  resonator  must  be  oriented  in 
such  a  way  that  the  precessing  magnetization,  M,  is  perpendicular  to  the  plane  of  the 
conductor  in  order  to  generate  the  maximum  signal.  The  main  characteristics  of  a 
resonator  are  the  frequency,  co,  at  which  it  resonates,  and  the  quality  factor,  Q,  which 
is  a  measure  of  the  circuit’s  energy  storage  efficiency.  At  resonance,  the  maximum 
peak  current  occurs  in  the  resonator  and  energy  is  alternatively  stored  as  magnetic 
energy  in  inductor,  and  then  as  electrostatic  energy  in  the  capacitor.  All  resonators, 
regardless  of  type,  must  be  tuned  to  the  resonant  frequency  of  the  nucleus  to  be 
studied.  They  also  have  to  be  matched  to  the  input  impedance  of  the  preamplifier  so 
that  the  maximum  power  is  transferred  from  the  resonator  to  the  receiver.  Tuning,  in 
general,  is  performed  via  a  tuning  circuit  made  up  of  a  fixed  and  variable  capacitor, 
while  impedance  matching  is  carried  out  using  a  variable  capacitor  only.  The  quality 
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factor  Q,  describing  the  capacity  of  the  resonator  to  store  magnetic  energy  without 
resistance  losses,  is  defined  as 

Q  =  maximum  energy  stored/average  energy  dissipated  per  radian  (1) 

or 

Q  -  (o/Afi)  (2) 


where  Aco  is  defined  as  the  full  3  dB  bandwidth.  The  definition  of  Q  given  by  (2)  lends 
itself  to  laboratory  measurements,  because  measuring  both  resonant  frequency  ca  and 
bandwidth  Aco  is  possible. 

A  good  RF  resonator  should  1)  resonate  at  the  desired  operating  frequency,  2)  have 
a  high  Q  factor,  3)  have  a  high  signal-to-noise  ratio  (SNR),  4)  have  a  good  filling 
factor,  5)  have  a  highly  homogeneous  Bt  field,  6)  have  minimum  resistance  losses,  7) 
be  able  to  withstand  the  applied  voltages,  8)  produce  minimal  electric  fields  in  the 
sample,  9)  have  minimum  interaction  with  the  rest  of  the  system,  10)  permit 
quadrature  excitation  and  reception,  and  11)  be  large  enough  to  accommodate  the 
imaging  sample.  This  project  demanded  a  new  efficient  resonator.  To  address  these 
needs,  I  designed,  built,  and  tested  two  RF  resonators,  one  being  a  double-tuned  whole 
body  ('H  -'^F)  birdcage  resonator  and  another  one  a  slotted  tube  resonator. 


2.2  Double-Tuned  Whole  Body  Birdcage  Resonator 

This  project  involved  '^F  and  'H  NMR  imaging  with  living  rats  in  a  single 
experiment,  so  it  was  highly  desirable  to  be  able  to  change  the  resonant  frequency 
without  retuning  the  resonator  or  disturbing  the  animal.  Thus,  1  chose  to  use  a  highpass 
birdcage  resonator,  since  it  provides  not  only  two  or  more  high  resonant  frequencies 
simultaneously,  but  also  the  greatest  possible  signal-to-noise  ratio  and  Bi  field 
homogeneity  [1].  A  birdcage  resonator  consists  of  a  set  of  N  (2,  4,  8,  16,  etc.)  copper 
wires  or  sheets  arranged  axially  on  the  surface  of  a  plastic  cylinder  and  connected  by 
high  quality  RF  capacitors  at  each  end.  In  this  configuration,  the  effective  current 
density  in  the  wires  or  sheets  varies  in  proportion  to  cos{(^,  where  ^  is  the  azimuthal 
angle  in  cylindrical  coordinates  [2],  and  there  exist  several  possible  modes  of 
resonance,  two  of  which  have  the  desired  sinusoidal  dependence  of  current  on  ^  [3]. 

2.2.1  Design  Theory 

Assuming  perfectly  conducting  copper  sheets.  Figure  1  shows  part  of  the 
simplified  circuit  model  of  a  birdcage  resonator  with  16  identical  legs  (copper  sheets), 
where  L,  is  the  self  inductance  of  the  /th  copper  sheet  and  C,  is  the  capacitor  connected 
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between  the  /th  and  the  (z+l)th  copper  sheet.  The  self  inductance  of  the  short  copper 
sheet  used  to  connect  the  capacitors  at  the  ends  is  neglected.  This  is  based  on  the  fact 
that  the  electric  current  in  the  resonator  is  longitudinal  (along  the  copper  sheet)  and 
maximal  at  the  ends  of  the  resonator. 


C,2  C..,  C,  C,.,  Ci.2 


Ci.2  Ci-i  Cj  Cj+1  Ci+2 


Figure  1  Part  of  the  equivalent  circuit  for  the  high-pass  double-tuned  whole  body 
birdcage  resonator. 


By  applying  Kirchhoff  s  Law  to  the  /th  loop  made  up  of  /th  and  (/+l)th  copper 
sheets  and  /th  capacitor  Cj  and  using  the  complex  number  notation  for  AC  current  and 
impedance,  we  can  get  a  linear  system  of  equations: 

=  0  (1) 

k  k  ^  i 

where  is  the  mutual  inductance  between  the  /th  and  kth  copper  sheets,  J,  is  the 

current  in  the  /th  capacitor,  is  the  current  in  the  Ath  copper  sheet,  /  is  the  imaginary 

number,  and  the  summation  is  carried  out  over  all  possible  values  of  index  k  from  0  to 
N-  1. 


There  are  many  ways  of  solving  Equation  (1),  either  analytically  or 
numerically.  The  solution  will  represent  the  various  resonant  modes  of  operation  of  the 
resonator.  This  circuit  model  is  a  simplified  one  and  does  not  take  into  account  the 
complications  that  exist  in  any  real  implementation,  such  as  the  self-inductance  of  the 
capacitors  and  the  inductance  of  the  end  rings.  Nevertheless,  the  general  form  of  the 
model  is  not  altered  by  these  factors,  and  the  qualitative  description  is  still  valid. 
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2.2.2  Computer  Simulation 

Instead  of  solving  the  above-described  complex  equation,  I  simulated  this 
circuit  model  on  a  PC  computer  using  P-Spice  [4],  The  complete  simulation  circuit 
model  is  shown  in  Figure  2.  The  desired  resonant  frequencies  are  188.22  MHz  for 
and  200. 16  MHz  for  'H  in  a  4.7  T  magnet.  The  inductance  produced  by  the  end  rings 
was  neglected  in  this  simulation  circuit,  since  it  only  accounts  for  ~  5%  of  the  total 
inductance.  However,  the  circuit  model  is  adequate  for  analyzing  the  essential 
frequency  characteristics  of  the  resonator. 


Figure  2  The  P-Spice  simulation  circuit  of  the  double-tuned  whole  body  birdcage  resonator. 

The  first  step  in  simulating  this  problem  was  to  estimate  the  inductance  of  the 
copper  sheets,  which  is  given  by  [3]. 

L  =  2/j[  In  {2hlb)  +  1/2  ]  (2) 
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where  L  is  the  mductance  in  nanoHenries,  h  is  the  height  in  centimeters,  and  b  is  the 
width  in  centimeters.  The  copper  sheets  have  a  height  of  16.5  cm  and  a  width  of  1.0 
cm,  with  a  thickness  of  ~  0.2  mm.  The  calculated  value  of  the  inductance  based  on 
Equation  (2)  was  113  nH. 

The  second  step  was  to  substitute  this  value  into  the  simulation  circuit,  apply  an 
AC  current  source  to  the  circuit,  vary  the  value  of  each  of  the  32  capacitors,  and 
observe  the  amplitudes  of  output  voltages  and  resonant  frequencies  of  the  circuit  until 
two  sharp  resonant  peaks  (188.22  MHz  for  '^F  and  200.16  MHz  for  ’H  in  a  4.7  T 
magnet)  with  the  maximum  output  voltages  were  achieved. 

The  capacitance  values  for  the  end  ring  capacitors  obtained  from  the  P-Spice 
simulation  study  are  tabulated  in  Table  1.  Theoretically  speaking,  all  the  capacitance 
values  should  be  equal  for  a  single-tuned  resonator  and  close  to  each  other  for  a 
double-tuned  resonator  in  order  to  generate  a  homogeneous  Bi  field.  The  circuit  model 
presented  here  is  a  simplified  one.  Therefore,  the  capacitance  values  obtained  are  not 
quite  close  to  each  other.  The  simulation  study  indicates  that  CIO,  C26,  Cll,  C27, 
C12,  and  C28  are  relatively  insensitive  to  '®F  resonant  frequency,  thus,  they  could  be 
connected  to  the  'H  tuning  circuit.  I  chose  to  place  the  tuning  circuit  in  parallel  to  CIO 
since  this  placement  provides  both  the  best  tuning  for  'H  mode  and  also  stabilizes  the 
mode. 


Table  1  Capacitance  values  determined  by  P-Spice  simulation  study  (unit:  pF) 


Fine  tuning  of  ’®F  mode  is  achieved  by  connecting  another  capacitor  in  parallel 
with  C6  as  shown  in  Figure  2.  When  '^F  tuning  capacitor  CF19T  varies  from  0  pF  to 
2000  pF,  '®F  resonant  peak  varies  from  (188.224  MHz,  14.727  V)  to  (184.871  MHz, 
16.606  V)  while  'H  resonant  peak  changes  only  from  (200.435  MHz,  13.602  V),  which 
is  a  clear  indication  that  each  mode  is  independently  tuned.  The  frequency  response 
characteristics  of  simulation  circuit  are  shown  in  Figure  3.  The  figure  shows  clean. 


10 


selective  excitations  of  the  individual  resonant  modes.  Simulation  assumes  a  high  Q 
circuit  (very  small  resistance). 


Figure  3  Frequency  response  characteristics  of  the  double-tuned  birdcage 
resonator  from  P-Spice  simulation  study 


2.2.3  Construction 

A  16-leg  high-pass  double-tuned  birdcage  resonator  was  constructed  according 
to  the  design  principles  and  the  simulation  results  described  above.  The  resonator  is 
tunable  to  both  188.2  MHz  and  200. 1  MHz  for  operation  at  ’"F  and  'H  frequencies  on  a 
4.7  T  magnet.  The  cylindrical  resonator  body  is  24.60  cm  in  length  and  has  an  outer 
diameter  of  24.00  cm.  The  outer  structure  of  the  resonator  is  made  of  a  5  mm  thick 
plexiglass  cylinder,  which  provides  mechanical  support  and  structural  stability  for  the 
resonator.  A  0.2  mm  thick  copper  sheet  shield,  coimected  to  the  ground,  is  mounted  on 
the  irmer  wall  of  the  outer  plexiglass  cylinder. 

The  critical  part  of  the  resonator  is  its  iimer  structure.  The  iimer  structure,  made 
from  a  5  mm  thick  plexiglass  cylinder,  has  an  outer  diameter  of  10.50  cm.  On  the  inner 
wall  of  the  plexiglass  cylinder  are  mounted  directly  16  copper  sheets  (legs)  with  1.5 
mm  gaps  and  32  capacitors.  The  size  of  the  umform  region  of  B\  magnetic  field 
generated  by  the  resonator  is  directly  related  to  the  number  of  legs.  The  more 
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numerous  the  legs,  the  more  uniform  the  Bx  magnetic  field  in  the  radial  direction.  I 
chose  to  build  a  birdcage  resonator  with  sixteen  legs  because  the  Bx  magnetic  field 
generated  should,  in  principle,  be  uniform  enough  for  my  applications  and  should 
produce  strong  NMR  signal.  The  legs  are  made  from  0. 1  mm  thick  copper  foil.  Each  of 
the  sixteen  legs  is  16.5  cm  in  length,  1.0  cm  in  width  in  the  central  part  and  1.5  cm  in 
width  at  the  two  ends  (Figure  4). 


- ■^2.0cm'r^ - 


I  I 


1.5  ciVi 

Figure  4  Geometry  of  the  legs  (copper  sheets). 


The  structure  of  the  two  end  rings,  i.e.,  the  two  circular  bands  at  the  ends  of  the 
resonator,  consisted  of  capacitors  and  was  made  identical  to  achieve  electrical  balance 
so  that  the  loading  effects  of  the  sample  could  be  reduced.  The  resonator  is  double- 
tuned,  and  is  driven  between  equivalent  points  on  the  two  rings  halfway  between  two 
adjacent  parallel  legs. 

The  resonator  was  tuned  by  individually  adjusting  the  CF19T  and  CHIT  at  each 
of  the  resonant  frequencies.  The  capacitively  coupled  tuning  loops  were  placed 
between  legs  11  and  13  and  legs  19  and  21  as  shown  in  Figure  2.  These  loops  were 
connected  by  standard  coaxial  cables  to  the  external  tuning  capacitors  CF19T  and 
CHIT  that  allowed  each  resonant  mode  to  be  inependently  tuned.  Special  care  was 
taken  to  ensure  that  there  was  little  tuning  interaction  between  the  two  resonant 
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frquencies  as  possible.  The  resonator  was  matched  capacitively  through  a  variahle 
capacitor  Cmatch  A  photo  of  the  resonator  is  shown  in  Figure  5. 


Figure  5  A  photo  of  the  16-leg  double-tuned  whole  body  birdcage  resonator. 


2.2.4  Bench  Testing 

The  laboratory  bench  testing  was  carried  out  using  a  network  analyzer 
(HP8752A).  To  test  the  resonator’s  performance,  it  was  essential  to  have  a  set  of 
standard  samples  that  simulated  the  type  of  in  vivo  samples  used  in  the  magnet.  Two 
NMR  phantoms  were  prepared  for  this  purpose,  which  roughly  had  the  same  loading 
effects  (lowering  Q  and  shifting  resonant  frequencies)  as  a  male  adult  rat  (~  250  g).  A 
444  ml  phantom  was  made  of  a  plastic  bottle  filled  with  a  solution  of  saline  and 
trifluoroacetic  acid  (TFA,  5%  v/v).  A  smaller  250  ml  phantom,  also  made  of  a  plastic 
bottle,  was  filled  with  a  solution  of  TFA  (5%  VA^)  and  CUSO4  (1  g/liter)  in  distilled 
water.  The  major  chemical  components  and  the  concentrations  of  the  phantoms  are 
listed  in  Table  2. 
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The  resonator’s  resonant  frequencies,  unloaded  and  loaded  Q  values  (quality 
factors)  and  Q  damping  for  each  of  the  resonant  modes  as  measured  by  the  network 
analyzer  are  presented  in  Table  3.  go  is  the  unloaded  Q  value  measured  with  the 
resonator  m  air  (without  phantom),  while  Qi^  is  the  Q  value  measured  when  the  saline 
phantom  was  placed  centrally  inside  the  resonator.  The  Q  damping  is  defined  as  the 
ratio  of  QolQh  ■  The  results  are  the  means  of  three  measurements.  The  loaded  Q  values 
are  affected  by  the  dimensions  of  the  phantom  and  its  position  in  the  resonator.  In 
general,  the  loaded  Q  values  are  dominated  by  the  effect  of  the  RF  lossy  sample.  The 
Q  values  decrease  as  the  dimensions  of  the  sample  increases. 


Table  2  Major  chemical  components  and  the  concentrations  of  the  444  ml 
saline  phantom. 


CHEMICAL  NAME 

CONCENTRATION 

(^1) 

NaCl 

6.90 

KCl 

0.375 

MgS04 

0.145 

CciCh 

0.185 

Na2HC03 

2.10 

Glucose 

1.80 

Table  3  Results  of  the  laboratory  bench  testing  of  the  birdcage  resonator, 
go  and  Ql  were  measured  with  the  444  ml  phantom. 


Resonant 

Mode 

Resonant 

Frequency 

(MHz) 

Unloaded 

Qo 

Loaded 

Ql 

Qo/Ql 

188.22 

14 

9 

1.56 

200.16 

13 

8 

1.63 

It  can  be  seen  from  Table  3  that  the  unloaded  go,  loaded  Qu  and  the  g  damping 
go  /gL  for  both  resonant  modes  are  similar.  The  relatively  low  g  values  may  be 
attributed  to  the  relatively  long  coaxial  cables  used  to  connect  the  tuning  and  match 
circuits.  Another  possible  reason  could  be  the  large  size  of  the  resonator  because  the 
longer  the  legs,  the  larger  their  resistance,  causing  increased  resistance  losses  and 
significantly  lowering  g  values. 
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One  of  the  consequences  of  the  low  Q  values  on  the  NMR  experiments  would 
be  to  increase  the  90°  pulse  width  for  a  given  RF  power  level.  This  can  be  verified  by 
the  following  equation  [4]: 


( 

PW  =  K 

V 


Vr 

Q 


y 


(3) 


where  PW  is  the  90°  pulse  width,  Vr  is  the  resonator  volume,  Q  is  the  quality  factor 
of  the  resonator,  and  the  factor  K  is  roughly  a  constant  for  all  resonator  designs.  This 
equation  iudicates  that  the  90°  pulse  width  is  linearly  proportional  to  the  square  root  of 
resonator’  volume  and  inversely  proportional  to  the  square  root  of  resonator’s  Q  value. 

2.2.5  Phantom  Imaging 

The  double-tuned  birdcage  resonator’s  NMR  performance  was  tested  on  the 
Omega  CSI  4.7  T  superconducting  magnet,  which  operates  at  a  proton  frequency  of 
200.106  MHz  and  a  fluorine  frequency  of  188.273  MHz.  The  phantom  was  centered  in 
the  resonator  and  the  resonator  was  positioned  in  the  isocenter  of  the  magnet,  with  its 
axis  being  aligned  with  the  direction  of  the  magnet’s  static  magnetic  field  Bo.  The  90° 
pulse  widths  were  determined  from  a  180°  null  of  Ihe  whole  sample.  Typical  90°  pulse 
widths  were  700  ps  for  proton  and  350  ps  for  fluorine.  Shimming  was  accomplished 
on  the  water  proton  FID  of  the  saline  solution  in  the  phantom,  with  the  resonator  tuned 
to  the  double  resonant  modes.  The  goal  was  to  obtain  a  Bo  of  homogeneous  strength 
over  the  entire  imaging  sample  volume.  Fourteen  shim  currents  were  adjusted  on  the 
proton  FID  throught  the  system  automated  shimming  utility  until  a  symmetrical  and 
long  FID  was  obtained.  The  proton  spectral  linewidth  at  half  height  after  shimming 
was  34  Hz,  which  was  sufficient  for  imaging  experiments. 

To  minimize  acquisition  time,  a  driven  equilibrium  spin-echo  pulse  sequence 
with  a  fairly  short  echo  time  (TE  =  8ms)  was  employed  so  that  the  T2  dephasing  effects 
were  not  significant.  Proton  transaxial  (transverse)  images  (Figure  6)  were  acquired  as 
a  3-D  data  set  with  a  repetition  time  TR  =  150  ms,  echo  time  TE  =  8  ms,  and  a  128  x 
64  X  8  matrix  size  over  100  mm  x  100  mm  x  200  mm  field  of  view  (FOV),  providing 
0.78  mm  x  1.56  mm  x  25  mm  digital  resolution.  To  improve  signal-to-noise  ratio 
(SNR),  two  acquisitions  were  averaged  for  each  image,  thus,  giving  a  total  acquisition 
time  of  2  min  33  sec.  Proton  coronal  (Figure  7(A)  and  (B))  and  sagittal  images  (Figure 
7  (C)  and  (D))  were  also  acquired  with  the  same  imaging  parameters  except  FOV  = 
200  X  200  mm. 

The  signal-to-noise  ratios  (SNRs)  for  the  images  with  2  averages  were 
measured  using  a  standard  routine.  The  values  for  signal  were  measured  from  a 
randomly  chosen  portion  of  images  by  a  square  crop  and  the  values  for  noise  were 
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taken  from  the  four  comers  of  the  image  background.  For  each  measurement,  the 
square  crop  covered  a  reasonably  large  area  of  the  image  or  background,  so  that  the 
data  obtained  were  unbiased.  The  results  of  SNR  measurements  for  the  proton  images 
are  listed  in  Table  4.  The  resonator’s  B\  field  homogeneity  was  estimated  by  visual 
inspection  of  the  images. 


Table  4  Results  of  SNR  meeisurements  for  the  proton  images. 


Upper- 

Left 

SNR 

Upper- 

^ght 

SNR 

Lower- 

Left 

SNR 

Lower- 

Right 

SNR 

Mean 

SNR 

SD 

Transaxial  Section 

88.0 

97.4 

93.0 

93.3 

92.9 

3.9 

Sagittal  Section 

279.8 

272.7 

278.7 

294.0 

281.3 

9.0 

Coronal  Section 

287.5 

114.3 

267.7 

125.8 

198.8 

91.4 

The  spectrometer  frequency  was  then  set  on  resonance  for  the  CF3  group  of 
TFA,  without  retiming  the  resonator.  The  90°  pulse  was  360  ms.  The  ’T  spin-echo 
images  of  the  444  ml  phantom  were  acquired  as  transaxial,  saggital  and  coronal 
projections  (Figure  8).  The  images  were  fluorine  density-weighted  with  TR  =  150  ms 
and  TE  =  8  ms  (using  driven  equilibrium  spin-echo  pulse  sequence)  and  were  averaged 
16  times  to  achieve  an  acceptable  SNR.  Thus,  each  image  was  acquired  in  2  min  33 
sec.  The  image  matrix  size  was  128  x  64  with  a  field  of  view  of  100  x  100  mm  for  the 
tranaxial  projection  and  200  x  200  mm  for  the  sagittal  and  coronal  projections, 
respectively. 


Table  5  Results  of  SNR  measurements  for  the  ’^F  images. 


Upper- 

Left 

SNR 

Upper- 

^ght 

SNR 

Lower- 

Left 

SNR 

Lower- 

Right 

SNR 

Mean 

SNR 

SD 

Transaxial  Section 

3.2 

3.3 

3.4 

3.5 

3.3 

0.1 

Sagittal  Section 

5.9 

6.7 

7.3 

6.9 

6.7 

0.6 

Coronal  Section 

4,3 

4.5 

3.8 

3.7 

4.1 

0.4 

The  image  SNRs  were  measured  in  the  same  way  as  the  proton  images  and  the 
results  are  presented  in  Table  5.  The  spin-lattice  relaxation  time  Ti  of  the  CF3  group 
was  estimated  by  applying  a  non-spatially  selective  inversion  recovery  (IR)  RF  pulse. 
Seven  different  IR  delays  (t)  increasing  in  the  range  between  800  ms  to  20  sec  were 
used  (Table  6).  The  Ti  value  was  calculated  based  on  the  Levenberg-Marquardt  three- 
parameter  fitting  algorithm  [5]  on  peak  intensity  values  (Table  6) . 
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Figure  6  Proton  transaxial  images  of  a  444  ml  phantom  acquired  using  a  4.7  T  magnet.  The 
images  shown  are  consecutive  slices.  Imaging  parameters:  TE  =  8  ms,  TR  =  150  ms,  slice 
thickness  =  2.5  cm,  NA  =  2,  FOV  =  100  x  100  mm  and  matrix  size  =  128  x  64. 
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Figure  7  Proton  coronal  (A  and  B)  and  sagittal  (C  and  D)  images  of  a  444  ml  phantom  acquired 
using  a  4.7  T  magnet.  Imaging  parameters:  TE  =  8  ms,  TR  =  150  ms,  slice  thickness  =  2.5  cm, 
NA  =  2,  FOV  =  200  x  200  mm  and  matrix  size  =  128  x  64. 
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Figure  8  MR  projection  images  of  a  444  ml  phantom.  Imaging  parameters:  TE  =  8  ms,  TR 
=  150  ms,  NA  =  16,  FOV  =  200  x  200  mm  (B  and  C)  ,  FOV  =  100  x  100  mm  (A)  and  matrix  size 
=  128  X  64..  A)  A  transaxial  projection  image.  B)  A  sagittal  projection  image.  C)  A  coronal 
projection  image. 
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Table  6  Results  of  T1  measurement  for  TFA  solution  and  delay  list  ( 7 )  of 
the  inversion  recovery  pulse  sequence. 


Function:  y  =  A*(l-(W+l)*exp(-T/Tl) 
A  =  1.92e+«6±26484 
W  =  0.026  ±0.042 
Tl=  2.3  ±0.19 
Standard  Error  =  41129.6 
Y  Standard  Deviation  =  31091.1 


Tau 

Intensity 

Calc 

Diff 

#sd 

20.000 

1923050.524 

1923610.434 

559.911 

0.018 

16.000 

1932423.278 

1921888.260 

-10535.018 

-0.339 

8.000 

1857983.273 

1859578.430 

1595.157 

0.051 

4.000 

1534324.010 

1567419.792 

33095.782 

1.064 

2.000 

1135584.863 

1085041.321 

-50543.542 

-1.626 

1.000 

587685.972 

637132.244 

49446.272 

1.590 

0.800 

545796.816 

522178.645 

-23618.170 

-0.760 

2.2.6  Discussion 

In  the  past,  the  complexity  of  birdcage  resonator  hindered  the  application  of 
traditional  equivalent  circuit  analysis  [3],  Here,  I  present  a  practical  approach  based  on 
computer  simulation  to  build  a  16-leg  double-tuned  whole  body  birdcage  resonator. 
The  resonator  has  two  resonant  frequencies  with  188.272  MHz  corresponding  to  ’®F 
and  200.106  MHz  to  'H.  Laboratory  bench  testing  of  the  resonator  demonstrated  that 
the  unloaded  Qs  for  both  proton  and  fluorine  resonant  modes  are  similar  and  relatively 
low.  The  loaded  Q  drops  by  about  36%  for  fluorine  and  39%  for  proton.  The  low  Q 
values  mean  inefficient  use  of  RF  power,  which  has  been  verified  by  the  NMR 
imaging  experiments.  As  a  result,  this  resonator  requires  the  use  of  a  high  RF  power 
level  that  could  lead  to  increased  heating  of  the  subject  under  investigation.  One  of  the 
possible  reasons  that  limits  the  Q  value  is  the  use  of  copper  shield  intended  to  shield 
the  sample  from  electric  field  interference.  In  fact,  the  shield  not  only  lowers  the  Q 
value  but  also  lengthens  the  pulse  width  [4]. 

The  Bi  field  was  found  to  be  somewhat  heterogeneous  across  usable  volume  of 
the  resonator  for  ’H  resonant  mode.  The  Bi  field  also  exhibited  nonuniformity  across 
the  phantom  for  'T  resonant  mode.  While  this  is  not  so  critical  for  the  sagittal  and 
transaxial  imaging  shown  here,  it  is  necessary  to  increase  the  RF  power  level  and 
sacrifice  some  pulse  width  to  achieve  an  acceptable  Bi  field  strength  for  the  coronal 
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imaging.  The  Bi  field  could  be  entirely  due  to  the  technical  difficulty  in  symmetrically 
positioning  and  adjusting  the  ring  capacitances  in  order  to  generate  an  ideal  sinusoidal 
distribution  of  current  around  the  resonator  at  the  higher  frequency.  Since  sensitivity 
of  the  resoantor  is  crucial  in  the  NMR  in  vivo  experiments,  1  believe  that  the 
dimensions  of  the  resonator  should  be  minimized  in  order  to  obtain  an  optimal  pulse 
width  for  the  particular  subject  of  investigation. 


2.3  Slotted  Tube  Resonator 

The  slotted-tube  resonator  was  originally  developed  to  efficiently  provide  'H 
decoupling  to  samples  in  high-field  magnets  [6].  The  common  slotted  tube  resonator 
consists  of  a  conducting  tube  with  two  slots  cut  symmetrically  along  both  sides  of  the 
tube.  The  surface  current  distribution  produced  by  the  slotted  tube  resonator  peaks 
near  the  edges  of  the  conducting  tube  adjacent  to  the  slots.  The  conducting  tube  of  the 
slotted  tube  resonator  is  opaque  to  RF  field  and  diverts  all  the  Bi  field  flux  through  the 
two  slots.  Thus,  the  aperture  angle  of  the  slots  must  be  optimized  in  order  for  the 
resonator  to  generate  a  homogeneous  field.  The  fundamental  feature  of  the  slotted  tube 
resonator  is  that  the  conducting  component  has  a  low  inductance,  and  hence,  large 
volume  resonators  can  be  constructed.  At  the  RF  frequencies  encountered  m  clinical 
MRI  systems,  this  means  that  the  slotted  tube  resonators  can  be  designed  to  contain 
samples  of  size  comparable  to  that  of  the  human  head.  The  major  applications  of  the 
slotted  tube  resonator  are  for  imaging  and  for  localized  spectroscopy  techniques 
demanding  a  uniform  flip  angle  over  the  entire  sample. 

2.3.1  Design  Theory 

The  circuit  model  is  illustrated  in  Figure  9.  Ls  and  R,  are  sample  inductance  and 
resistance,  respectively.  This  circuit  model  can  be  further  simplified  in  terms  of 
equivalent  circuit  (Figure  10)  with 


i(L,C,co^  -1) 
coiC,-C,(L,C,co^-l)) 


(4) 


where  co  is  the  resonant  frequency  of  the  circuit  and  i  is  the  imaginary  number.  Thus, 
the  tuning  impedance  Zj  is  given  by  the  following  expression: 


‘  "  C^o)  " 

.  Q  - 1)  - u>(C3  -  C, iL,C,co^  - 1)) 

ojC,{C, -C,{L,C,C^  -1)) 


(5) 
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Figure  9  The  circuit  model  of  the  switchable  slotted  tube  resonator.  R*  is  the  sample 
Resistance  and  the  sample  inductance.  R3  is  the  equivalent  resistance  of  the  tuning  loop. 


Figure  10  Equivalent  circuit  model  of  the  slotted  tube  resonator. 


When  the  resonator  operates  at  the  higher  frequency  ty  =  tyt,,  the  tuning 
impedance  Zj  is  a  fimction  of  the  total  capacitance  of  the  tuning  circuit  and  is  given  by 
the  expression 


(6) 
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where 


_C,{L,C,a>l-l)-C, 
'  L,C,col-l 


(7) 


Thus,  the  tuning  is  achieved  by  adjusting  all  the  capacitors. 

When  the  resonator  operates  at  the  lower  frequency  (0=  0)i ,,  \he,  series  Ls  C3 
circuit  has  a  zero  impedance  and  the  total  tuning  impedance  is  reduced  to 


CjU)i 


(8) 


Thus,  tuning  is  obtained  by  adjusting  Ci  only. 

2.3.2  Construction 

A  switchable  slotted  tube  resonator  was  constructed  based  on  the  above  design 
theory  using  two  copper  sheets,  which  cover  an  80°  arc  [9]  for  optimal  Bi  field 
homogeneity.  The  cylindrical  resonator  body  is  23.5  cm  long  with  outer  diameter  25.4 
cm.  The  outer  structure  of  the  resonator  is  made  of  a  5mm  thick  plexiglass  cylinder. 
The  inner  structure  is  made  of  a  3  mm  thick  plexiglass  cylinder,  having  a  small  outer 
diameter  of  7.6  cm  to  increase  the  filling  factor  and  reduce  the  pulse  width.  The  two 
ends  of  the  resonator  are  covered  by  two  6  mm  thick  plexiglass  boards  stabilized  by  6 
nonmagnetic  screws  on  each  end,  respectively.  The  two  boards,  along  with  the  outer 
structure,  provide  very  stable  mechanical  support  for  the  resonator.  Three  capacitors  of 
20.0  pF  were  positioned  symmetrically  between  the  guard  rings  and  the  copper  sheets 
to  obtain  uniform  surface  current  distribution. 

The  resonator  has  a  usable  length  of  15.5  cm,  sufficient  for  in  vivo  adult  rat 
NMR  experiments.  The  two  guard  rings  and  two  copper  sheets  are  made  from  0.3  mm 
thick  copper  foil.  The  width  of  each  of  the  guard  rings  is  1.0  cm  and  there  is  a  spacing 
of  approximately  0.5  cm  between  the  guard  ring  and  the  vertical  copper  sheets,  which 
are  joined  with  capacitors.  The  copper  sheets  have  a  length  of  12.5  cm  and  their  width 
covers  about  80°.  The  resonator  is  tunable  to  '^F  (188.22  MHz)  and  ‘H  (200.16  MHz), 
and  is  driven  between  one  copper  sheet  and  one  guard  ring. 

The  tuning  of  the  resonator  was  achieved  by  three  variable  capacitors  Ci,  C2 
and  C3,  which  are  V2145H  family  high  Q  RF  capacitors  (Voltronics  Corporation, 


23 


Denville,  New  Jersey).  Ci  and  C2  (NMTM38GE)  have  a  capacitance  range  of  1.0  ~ 
38.0  pF,  C3  (NMTM120CE)  has  a  capacitance  range  of  2.0  ~  120.0  pF,  and  Cm 
(NMQM22GE)  has  a  capacitance  range  of  1.0  ~  22.0  pF.  '®F  resonance  was  obtained 
by  adjusting  Ci  and  'H  resonance  was  achieved  by  adjusting  Ci,  C2,  and  C3, 
respectively.  The  resonator  was  matched  capacitively  through  Cm-  A  photo  of  the 
resonator  is  shown  in  Figure  11. 


Figure  10  A  photo  of  the  slotted  tube  whole  body  resonator. 


2.3.3  Bench  Testing 

The  resonator  was  tuned  and  matched  on  the  laboratory  bench  with  the  help  of  a 
network  analyzer  (HP8752A).  The  unloaded  and  loaded  Q  values  were  measured  with 
the  aid  of  the  network  analyzer  from  the  full  width  at  half  height  of  the  resonance 
when  the  resonator  was  perfectly  tuned  and  matched. 

The  resonant  frequencies,  unloaded  Qa  and  loaded  0l  values  and  Q  damping  for 
each  of  the  resonant  modes  are  presented  in  Table  7.  Qq  is  the  unloaded  Q  value 
measured  with  the  resonator  in  air  (without  phantom),  while  0l  is  the  Q  value 
measured  when  the  phantom  was  placed  centrally  inside  the  resonator.  The  Q  damping 
is  defined  as  the  ratio  of  Qo  /Qi^ .  The  results  are  Ihe  means  of  three  measurements. 
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Table  7  Results  of  the  laboratory  bench  testing  of  the  slotted  tube  resonator. 
Qo  and  Ql  were  measured  with  the  250  ml  phantom. 


Resonant 

Mode 

Resonant 

Frequency 

(MHz) 

Unloaded 

Qo 

Loaded 

Ql 

Qo/Ql 

'"F 

188.22 

97.5 

54.2 

1.80 

mm 

200.16 

63.8 

53.6 

1.19 

As  can  be  seen  from  Table  7,  the  unloaded  Qo  for  '®F  is  much  higher  than  that 
for  'H  1.5  times),  and  the  loaded  Ql  for  both  resonant  frequencies  is  essentially  same. 
The  Q  damping,  Qo  /Ql,  for  '®F  is  1.5  times  higher  than  that  for  ’H.  This  means  that 
the  RF  lossy  phantom  has  a  much  stronger  effect  on  the  'T  resonance  than  on  the  ’H 
resonance.  The  high  Q  values  could  be  the  result  of  several  factors,  including  a 
relatively  small  usable  volume,  minimized  length  of  connecting  wires,  and  a  structure 
without  a  copper  shield.  The  liigh  Q  values  guarantee  that  the  resonator  has  a  higher 
excitation  efficiency,  i.e.,  the  90°  pulse  width  for  a  given  RF  power  level  is  short  or  the 
magnitude  of  the  field  generated  with  a  given  RF  power  level  is  high.  The  results  of 
the  laboratory  bench  testing  also  indicated  that  the  resonator  was  easy  to  tune  and 
match  and  both  resonant  frequencies  were  stable  and  immune  to  external 
electromagnetic  interference. 


2.3.4  Phantom  Imaging 

The  NMR  performance  of  the  slotted  tube  resonator  was  tested  on  the  Omega 
CSI  4.7  T  superconducting  horizontal  bore  magnet  operating  at  a  proton  resonant 
frequency  of  200.106  MHz  and  a  fluorine  resonant  frequency  of  188.273  MHz.  The  in 
vitro  resonator  sensitivity  for  proton  and  fluorine  was  determined  using  the  250  ml 
phantom,  which  was  filled  with  a  solution  of  TFA  (5%v/v)  and  CUSO4  (Ig/liter)  in 
distilled  water. 

Measurements  of  90°  pulse  width  were  performed  from  a  180°  null  of  the  whole 
sample.  Typical  90°  pulse  widths  were  around  170  |j,s  for  proton  and  110  jxs  for 
fluoiine.  Shimming  was  performed  on  the  water  proton  FID  of  the  saline  solution  m 
the  phantom  with  resonator  tuned  to  the  'H  resonance,  to  a  typical  spectral  linewidth  of 
32  Hz.  The  driven  equilibrium  spin-echo  pulse  sequence  was  used  for  the  imaging 
experiment.  Imaging  parameters  were:  repetition  time  TR=  150ms,  echo  time  TE  = 
8ms,  with  a  data  set  of  128  x  64  x  8  with  128  real  and  128  imaginary  data  points,  64 
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phase  encoding  steps,  and  8  slices.  The  field  of  view  (FOV)  was  80  mm  x  80  mm  for 
transaxial  images,  100  mm  x  100  mm  for  coronal,  and  200  mm  x  200  mm  for  sagittal 
images.  The  images  were  acquired  with  one  excitation,  thus  giving  a  total  acquisition 
time  of  1.16  min  (Figures  12,  13,  14,  and  15).  The  image  signal-to-noise  ratios  (SNRs) 
were  measured  as  described  before.  The  results  of  SNR  measurements  for  the  proton 
images  are  listed  in  Table8. 


Table  8  Results  of  SNR  measurements  for  the  proton  images. 


Upper- 

Left 

SNR 

Upper- 

Right 

SNR 

Lower- 

Left 

SNR 

Lower- 

Right 

SNR 

imiiiii 

SD 

Transaxial  Section 

389.3 

199.5 

897.8 

666.3 

177.1 

Sagittal  Section 

494.2 

517.5 

478.8 

485.3 

28.1 

Coronal  Section 

190.3 

30.2 

164.6 

102.0 

87.9 

Following  the  'H  imaging  experiment,  the  resonator  was  retuned  m  place  to 
188.27  MHz  and  corresponding  '^F  imaging  experiment  was  carried  out.  The  90°  pulse 
was  108  ps  and  the  spectral  linewidth  after  shimming  was  46  Hz.  The  spin-echo  (SE) 
images  were  acquired  with  the  imaging  parameters:  TR  =  150  ms,  TE  =  8  ms,  FOV  = 
100  X  100  mm,  NA  =  16,  and  matrix  size  =  64  x  32.  Figure  16  shows  four  fluorine  SE 
transaxial  images.  Figure  17  (A)  shows  a  proton  transaxial  EPl  image  and  Figure  17 
(B)  shows  the  corresponding  fluorine  EPI  image  of  the  same  phantom.  The  image 
SNR  measurements  were  performed  in  the  same  way  as  the  proton  images  and  the 
results  are  given  in  Table  9. 


Table  9  Results  of  SNR  measurements  for  the  fluorine  imasges. 


Upper- 

Left 

SNR 

Upper- 

Right 

SNR 

Lower- 

Left 

SNR 

Lower- 

Right 

SNR 

nmn 

SD 

Transaxial  Section 

3.1 

2.9 

2.6 

3.0 

0.2 

Sagittal  Section 

4.2 

5.1 

4.0 

4.2 

4.4 

0.5 

Coronal  Section 

4.2 

4.3 

4.7 

4.6 

4.5 

0.2 

26 


Figure  12  Proton  transaxial  images  of  a  250  ml  phantom  acquired  using  a  4.7  T  magnet.  The 
images  shown  are  the  first  4  slices.  Imaging  parameters:  TE  =  8  ms,  TR  =  150  ms,  slice  thickness 
=  12.5  mm,  NA  =  1,  FOV  =  80  x  80  mm  and  matrix  size  =  128  x  64. 


Figure  13  Proton  transaxial  images  of  a  250  ml  phantom  acquired  using  a  4.7  T  magnet.  The 
images  shown  are  the  last  4  slices.  Imaging  parameters:  TE  =  8  ms,  TR  =  150  ms,  sUce  thickness  = 
12.5  mm,  NA  =  1,  FOV  =  80  x  80  mm  and  matrix  size  =  128  x  64. 
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Figure  14  Proton  coronal  images  of  a  250  ml  phantom  acquired  using  a  4.7  T  magnet.  The  images 
shown  are  the  middle  4  slices.  Imaging  parameters:  TE  =  8  ms,  TR  =  150  ms,  slice  thickness  =  10 
mm,  NA  =  1,  FOV  =  100  x  100  mm  and  matrix  size  =  128  x  64 
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Figure  15  Proton  sagittal  images  of  a  250  ml  phantom  acquired  using  a  4.7  T  magnet  The 
images  shown  are  the  middle  4  slices.  Imaging  parameters:  TE  =  8  ms,  TR  =  150  ms,  slice 
thickness  =  10  mm,  NA  =  1,  FOV  =200  x  200  mm  and  matrix  size  =  128  x  64 
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Figure  16  Fluorine  transaxial  images  of  a  250  ml  phantom  acquired  using  a  4.7  T  magnet. 
The  images  shown  are  the  middle  4  slices.  Imaging  parameters:  TE  =  8  ms,  TR  =  150  ms, 
slice  thickness  =  10  mm,  NA  =  16,  FOV  =  80  x  80  mm  and  matrix  size  =  64  x  32 
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Figure  17  A)  A  proton  transaxial  EPI  image  of  a  250  ml  phantom  acquired  using  a  4.7  T 
magnet.  B)  A  fluorine  transaxial  EPI  image  of  the  same  phantom.  Imaging  parameters:  TE 
=  8  ms,  TR  =  1500  ms,  slice  thickness  =  10  mm,  NA  =  1  CH).  NA  =  4  (^^F),  FOV  =  80  x 
80  mm  and  matrix  size  =  32  x  32 


2.3.5  Discussion 

A  theoretical  approach  to  design  a  swithable  slotted-tube  resonator  was 
presented.  The  resonator  has  the  capability  of  being  continuously  tuned  in  the 
frequency  range  of  150  MHz  ~  220  MHz.  At  each  particular  frequency,  a  sharp 
resonance  can  be  easily  obtained  with  tuning  and  matching  circuits.  The  resonator  has 
high  Q  values  and  short  90 “pulse  widths  for  both  proton  and  fluorine,  which  results  in 
improved  sensitivity  and  makes  in  vivo  tumor  p02  studies  possible. 

A  comparison  between  Tables  4  and  8,  and  Tables  5  and  9  indicates  that  the  'H 
SNRs  for  the  slotted  tube  resonator  are  about  3  times  higher  than  those  for  the 
birdcage  resonator  and  the  SNRs  are  virtually  same  for  the  both  resonators. 
Several  possible  reasons  could  cause  this  discrepancy.  First  of  all,  the  slotted  tube 
resonator  has  a  higher  filling  factor  and  less  resistance  loss  than  the  birdcage  resonator, 
and  these  would  enhance  the  RF  transmitting  efficiency.  Secondly,  the  slotted  tube 
resonator  used  the  small  phantom  for  the  SNR  measurements  and  hence,  the  amount  of 
energy  received  per  nucleus  was  higher  than  that  for  the  large  phantom  for  a  given  RF 
power.  These  two  reasons  could  result  in  higher  'H  SNRs  for  the  slotted  tube 
resonator.  As  to  the  relatively  low  '®F  SNRs  for  the  slotted  tube  resonator,  this  could 
be  entirely  due  to  the  low  fluorine  nucleus  concentration  in  the  small  phantom 
solution.  Another  thing  worth  mentioning  here  is  that  when  acquiring  EPI  images,  the 
1  kW  amplifier  had  to  be  tuned  on  in  order  to  get  strong  '®F  signals  for  both  resonators. 

The  Bi  field  homogeneity  was  found  to  be  excellent  for  proton  resonance  and 
reasonably  good  for  fluorine  resonance  across  the  usable  volume  of  the  resonator.  The 
Bt  field  homogeneity  was  sufficient  for  this  project. 


3.  Software  Development 


3.1  Program  to  Assess  the  Goodness  of  T1  Relaxation  Curve  Fit 

EPI  is  used  as  a  basic  building  block  and  combined  with  additional  spin 
preparation  sequences.  To  measure  spin-lattice  relaxation  time  (Tl)  and  hence  /’Oj 
map,  a  pulse  burst  saturating  (PBSR)  pulse  train  is  applied  immediately  prior  to  the 
EPI  sequence.  This  leads  to  some  loss  of  dynamic  range  in  the  results,  but  substantially 
shortens  the  time  required  for  a  Tl  experiment  and  thus,  improves  the  temporal 
resolution.  Fourteen  data  sets  are  acquired  with  delay  time  in  a  geometric  progression 
ranging  from  200  ms  to  90  sec.  Since  the  longest  and  shortest  delays  are  alternated 
during  data  acquisitions  to  reduce  systematic  bias,  it  is  necessary  to  restore  the  images 

to  the  correct  order,  i.e.,  from  the  shortest  to  the  longest  delay  (200  ms . ^90  sec) 

prior  to  the  curve  fit 

In  order  to  obtain  the  spatial  distribution  (map)  of  p02  in  breast  tumor,  we  have 
to  solve  for  the  R1  map  of  injected  HFB  first.  This  is  done  by  fitting  each  voxel  of  the 
fourteen  images  to  a  three  parameter  exponential  relaxation  model: 


yn  0^  J)  =  J)  •  [1  -  (1  +  ^)  •  exp(-Rl(/,  j)  ■  t„ )] 

(«  =  1,2, . ,14)  (9) 

(/,;  =  i,2, . 32) 


by  the  Levenberg-Marquardt  least-squares  algorithm,  where  is  the  measured 

signal  intensity  corresponding  to  delay  time  T„{n\h  image)  for  voxel  is 

the  fully  relaxed  signal  intensity  amplitude  of  voxel  IF  is  a  dimensionless 

scaling  factor  allowing  for  imperfect  signal  conversion,  and  R\{i,j)  is  the  relaxation 
rate  of  voxel  {i,j)  in  unit  of  sec“^ . 

Good  Tl  relaxation  curve  fit  is  crucial  in  '®F  PBSR-EPI  oximetry.  It  is 
relatively  easy  to  obtain  excellent  curve  fits,  provided  SNR  >  10  for  most  intense 
signals.  To  assess  the  goodness  of  curve  fit,  I  wrote  a  Unix  based  NMR-shell  script 
program  to  display  the  Tl  relaxation  curve  of  individual  voxels  of  PBSR-EPI  images. 
The  graphical  representation  of  the  data,  relaxation  model,  and  goodness  of  fit 
provides  a  quick  and  convenient  way  of  assessing  the  quality  of  the  data.  Based  on  the 
rms  error  of  fit,  the  goodness  of  curve  fit  can  be  classified  into  three  categories:  good, 
intermediate,  and  poor. 
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For  typical  good  relaxation  data,  the  global  least  squares  minimization  is 
achieved  with  a  very  small  rms  error.  These  data,  related  to  strong  ’®F  signal,  represent 
those  voxels  within  the  region  of  interest  (ROI).  For  typical  intermediate  data,  the 
global  least  squares  minimization  is  achieved  with  a  relatively  large  rms  error  due  to  a 
signal  bad  data  point,  which  could  lead  to  a  big  T1  error.  To  improve  the  relaxation 
cmve  fit,  a  threshold  is  applied  to  the  raw  data,  which  eliminates  those  data  points 
>  3<t  from  the  curve.  The  remaining  data,  thus,  provide  a  better  curve  fit.  In  applying 
the  threshold,  only  one  data  point  is  eliminated  from  the  raw  data  per  curve.  These 
relaxation  data  represent  those  voxels  within  the  ROI,  but  have  relatively  strong  ’®F 
signals.  In  general,  they  give  good  T1  values.  For  poor  relaxation  data,  the  global  least 
squares  minimization  either  fails  or  converges  slowly  with  a  very  large  rms  error. 
These  relaxation  data,  associated  with  those  voxels  outside  the  ROI,  represent  either 
very  week  ’®F  signals  or  background  noise.  They  give  umeliable  T1  values. 


3.2  Program  to  Compute  HFB  Clearance  Rate  Map 

In  contrast  to  traditional  PFCs,  which  exhibit  excessive  tissue  retention,  HFB 
clears  from  tissue  relatively  rapidly.  Since  HFB  is  a  non-ionic  freely  diffusable  tracer, 
clearance  will  provide  an  indication  of  relative  breast  tumor  blood  flow.  To  compute 
the  HFB  clearance  rate  map,  I  use  a  single  parameter  exponential  decay  model: 


aj(O  =  0,,/(O)exp(-?/7;.^) 


(10) 


where  Qi  j  {t)  is  the  signal  intensity  for  voxel  {i,  j)  at  time  t,  0,  y(O)  is  the  signal 
intensity  for  voxel  (/,  j)  at  time  zero,  and  T^  j  is  the  clearance  rate  of  voxel  (/,  j).  In 
practice,  however.  Equation  (10)  is  transformed  into  a  logarithmic  form: 

ln(C,,(0)=ln(a,,(0))-?/7;.,  (11) 


Based  on  above  theory,  I  wrote  a  Unix  based  NMR-shell  script  program,  which 
computes  the  clearance  rate  T^j  of  individual  voxels  by  fitting  the  signal  intensity  of 

PBSR-EPI  images  to  Equation  (11).  A  clearance  rate  map  is  then  produced  by 
displaying  T^j  as  a  color-coded  image.  This  program  is  capable  of  determining  the 

global  and  regional  clearance  rate  of  HFB.  Source  code  of  the  program  is  provided  in 
Appendix. 
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4.  Conclusions 


The  construction  of  high  quality  RF  resonators  and  development  of  two 
important  programs  have  established  a  solid  foundation  for  this  project  and  made 
monitoring  global  and  regional  breast  tumor  oxygenation  using  'T  MRI  of  HFB 
possible.  PBSR-EPl  oximetry  of  HFB  will  provide  important  insights  into  breast  tumor 
physiology  and  therapy.  I  am  currently  analyzing  tumor  oxygenation  data  using  these 
two  programs  and  the  results  will  be  available  soon. 
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6.  Appendix 


6.1  Source  code  to  Compute  HFB  Clearance  Rate  Map 

#  Omega  Sourceable  Script 

#  File:  clearance_rate 

#  Subprograms  called:  clearance_rate_subl  and  clearance_rate_sub2 

#  Author:  Yulin  Song 

#  Status :  OK 

#  Time:  June  1998 


onintr  quit 
clear 
bigwin 
clear 


@  choice  =  1 
while  ( ($choice 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo 
echo  "" 
echo  " 
echo  "" 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  "" 


1)  I  I  ($choice  ==  2)11  ($choice  ==  3)  |  |  {$choice  = 


CLEARANCE  RATE  COMPUTATION  PROGRAM 
MAIN  MENU 

1.  Purposes  of  the  Program 

2.  Usage  of  the  Program 

3.  Display  T  Fitting  Curve 

4 .  Generate  T  Maps 

5.  Exit 
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II II 


if  (  $choice  ==  1  )  then 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 

echo  "  PURPOSES  OF  THE  PROGRAM 

echo  "" 

echo  "  ================================:=========:: 


echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
Echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 
echo  " 


1)  To  calculate  the  HFB  clearance  rate  {  biological 
half-life  )  T  of  each  individual  voxels  based  on 
mathematical  model: 

q{t)  =  q{0) *exp (-t/T) 

where  q(t)  is  the  signal  intensity  at  time  t,  q(0) 
the  signal  intensity  at  time  0 

(immediately  after  IT  injection  of  HFB  ),  and  T  the 
clearance  rate  in  min. 


2)  To  display  the  T  fitting  curve  of  each  individual 
voxels  based  on  the  logarithmic  transformation  of 
the  above  model: 

In  {q{t) )  =  In  (q{o) )  -  t/T 


3)  To  generate  the  T  maps  {  clearance  rate  maps  ) 
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echo  "  I 

echo  "  ==========  =======zz==========z=====  =  =  z==  =  =  : 

echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 

echo  -n  "  Please  press  RRTURN  to  continue  ->  : 

set  blank  =  $< 

clear 

else  if  {  $choice  ==  2  )  then 
clear 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 

echo  "  USAGE  OF  THE  PROGRAM 

echo  "  - - - - - - - - — 

echo  "  I 

echo  "  I 

echo  "I  1.  Run  ' iraqautosong '  program  one  hour 
echo  "  I  injection  of  HFB. 

echo  "  I 

echo  "I  2.  Move  the  data  with  the  longest  delay 
echo  "  I  program's  directory, 

echo  "  I 

echo  "1  3.  Rename  the  data  as  ' clearance_raw. 1 ' 

echo  "  I  you  prefer,  but  it  must  end  with  ext 

echo  "  I 

echo  "I  4.  Rerun  'iraqautosong'  every  30  min  or 
echo  "  I  time  interval  you  prefer, 

echo  "  I 

echo  "I  5.  Repeat  steps  2,  3,  and  4.  But  rename 
echo  "  I  end  with  '.2',  '.3',  '.4',  etc. 

echo  "  I  until  11  data  sets  have  been  acquire 

echo  "  I  renamed, 

echo  "  I 

echo  "I  6.  Run  this  program  and  follow  the  inst 
echo  "  I  the  screen, 

echo  "  I 

echo  "  I 

echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 


Run  'iraqautosong'  program  one  hour  after  IT 
injection  of  HFB. 

Move  the  data  with  the  longest  delay  to  this 
program's  directory. 

Rename  the  data  as  ' clearance_raw. 1 '  or  any  name 
you  prefer,  but  it  must  end  with  extension  '.1'. 

Rerun  'iraqautosong'  every  30  min  or  any  reasonable 
time  interval  you  prefer. 

Repeat  steps  2,  3,  and  4.  But  renamed  data  must 
end  with  '.2',  '.3',  '.4',  etc. 

until  11  data  sets  have  been  acquired  and  properly 
renamed. 

Run  this  program  and  follow  the  instructions  on 
the  screen. 
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It  II 


echo 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 

echo  -n  "  Please  press  RETURN  to  continue  ->  : 

set  blank  =  $< 

clear 

else  if  (  $choice  ==  3  )  then 
clear_rate_subl 
echo  "" 

echo  -n  "  Please  press  RETURN  to  continue  ->  : 

set  blank  =  $< 

clear 

else  if  {  $choice  ==  4  )  then 
clear 
sub2 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 

echo  -n  "  Please  press  RETURN  to  continue  ->  : 

set  blank  =  $< 
clear 
else 

smallwin 
echo  "" 


Please  press  RETURN  to  continue  -> 
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echo  "  Bye,  Bye  ! 
echo  "" 
echo  "" 
endif 


#  end  of  while 
end 


unset  blank 
quit ; 
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<  * 


,  i  ^  \ 


Omega  Sourceable  Script 
Filename:  clear_rate_subl 

This  is  the  first  subroutine  called  by  main  program  clear_rate 
Purpose:  To  display  the  T  (clearance  rate)  fitting  curve 
and  compute  the  clearance  rate 
Author:  Yulin  Song 
Time:  June  1998 
Status :  OK 

Potential  Problem:  The  model  involed  ln(x)  function  may  not  result  in  a 

convergent  fitting. 

In  this  case,  two  models  could  be  used: 

Y  =  SLOPE*x  +  INTERCEPT  or 

Y  =  Y0*exp (-x/T) 


set  file_dir  =  "$home/yulindata" 
set  fn  = 
set  n  = 

set  time  interval  = 


set  fn_temp  = 
set  n_temp  = 

set  time_interval_temp  = 
set  loop  =  N 

while  {  (  $loop  ==  N  )  II 

clear 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  " 
echo  "" 
echo  "" 


(  $loop  ==  n  )  ) 


echo  "" 
echo  "" 
echo  "" 


T  FITTING  CURVE  PARAMETER  MENU 


echo 

If 

- - - 

- ..... - 

— ... — ... — - - - - 

ff 

echo 

If  ff 

echo 

ff 

1. 

PROGRAM  NAME: 

=  Clearance  rate 

ft 

echo 

tf  ff 

echo 

ff 

2. 

CURRENT  DIRECTORY: 

=  $file  dir 

ff 

echo 

ff  ff 

echo 

ff 

3. 

FILENAME  BASE: 

=  $fn 

tf 

echo 

ff  ff 

echo 

ff 

4. 

NUMBER  OF  FILES: 

=  $n 

ft 

echo 

tf  ff 

echo 

ft 

5. 

TIME  INTERVAL  FILE 

=  $time  interval 

ft 

echo 

tf  ft 

echo 

ff 

— 

— 

- - 

ff 
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echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 


echo 

-n  " 

Please  enter  FILENAME  BASE 

->  " 

set  fn  = 

$< 

if  ( 

$fn 

j  = 

""  )  then 

set 

fn_ 

temp  =  $fn 

else 

set 

fn 

=  $fn  temp 

endif 

clear 

echo 

tf  tl 

echo 

II  II 

echo 

M  tl 

echo 

II  II 

echo 

II  II 

echo 

II  It 

echo 

It  ft 

echo 

It  It 

echo 

II  II 

echo 

II  II 

echo 

It 

T  FITTING  CURVE 

PARAMETER  MENU 

echo 

II  It 

echo 

echo 

echo 

It  It 

tl 

It  II 

- - - - - - — 

- - - - - 

echo 

It 

1.  PROGRAM  NAME: 

=  Clearance  rate 

echo 

M  II 

echo 

II 

2.  CURRENT  DIRECTORY: 

=  $file_dir 

echo 

II  It 

echo 

II 

3.  FILENAME  BASE: 

=  $fn 

echo 

II  II 

echo 

II 

4 .  NUMBER  OF  FILES : 

=  $n 

echo 

It  It 

echo 

It 

5.  TIME  INTERVAL  FILE 

=  $time  interval 

echo 

echo 

echo 

It  It 

It 

It  It 

echo 

II  II 

echo 

II  II 

echo 

II  II 

echo 

It  II 

echo 

II  II 

echo 

II  II 
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Please  enter  the  NUMBER  OF  FILES  -> 


echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  -n  " 
set  n  =  $< 
if  (  $n  !=  ""  )  then 
set  n_temp  =  $n 

else 


set  n  =  $n_temp 

endif 
clear 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  " 
echo  "" 
echo  "" 
echo  " 
echo  "" 


echo  "  1 

echo  "" 

echo  "  2 

echo  "" 

echo  "  3 

echo  "" 

echo  "  4 

echo  "" 

echo  "  5 

echo  "" 

echo  "  - 

echo  "" 


echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 


T  FITTING  CURVE  PARAMETER  MENU 


PROGRAM  NAME: 
CURRENT  DIRECTORY: 
FILENAME  BASE: 
NUMBER  OF  FILES: 
TIME  INTERVAL  FILE 


Clear ance_r ate 
$file_dir 
$  fn 
$n 

$time  interval 
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J  I 


«  ll  t  *- 


echo  "" 
echo  "" 
echo  "" 
echo  "" 

echo  -n  "  Please  enter  the  TIME  INTERVAL  FILE  -> 

set  time_interval  =  $< 

if  {  $time_interval  !=  ""  )  then 

set  time_interval_te;mp  =  $time_interval 

else 


set  time_interval  =  $time_interval_temp 


endif 
clear 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  " 
echo  "" 
echo  "" 
echo  " 
echo  "" 
echo  " 
echo  "" 
echo  " 
echo  "" 
echo  " 
echo  "" 
echo  " 
echo  "" 
echo  " 
echo  "" 
echo  " 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 
echo  "" 


T  FITTING  CURVE  PARAMETER  MENU 


1.  PROGRAM  NAME: 

2.  CURRENT  DIRECTORY: 

3.  FILENAME  BASE: 

4.  NUMBER  OF  FILES: 

5.  TIME  INTERVAL  FILE 


=  Clearance_rate 
=  $file_dir 
=  $fn 
=  $n 

=  $time  interval 


M 


ft 


If 


It 


It 


II 


11 
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echo  -n  "Are  the  entries  correct  (Y  for  correct  and  N  for  wrong)  ->  : 

set  loop  =  $< 

end 

#  end  of  while 


stackrm  * 
get  $fn.l 
dis 

@  fit_another_voxel  =  1 
while  (  $fit_another_voxel  ) 

#  the  first  while  statement, 
echo  "" 

echo  "  Please  enter  the  coordinates  of  the  voxel  to  be  plotted:  " 
echo  "" 

echo  -n  "  X  =  ->  " 
set  col  =  $< 

cal  -o  col  $col  1  sub  #  col  =  col  -  1 

echo  "" 

echo  -n  "  Y  =  ->  " 
set  row  =  $< 
echo  "" 

stackrm  * 
rm  -f  %tmp 


@  ex  =  $n 

while  (  $ex  >=  1  ) 

#  get  the  files  in  descs'.nting  order  i.e.  in  decreasing  signal 

#  intensity  order 

get  $fn.$ex 
@  ex — 

end 


#  cp  $fn.l  junk 

#  get  junk 

header  -i  dO  dimO  #  no  of  columns:  acquisition  direction 

header  -i  dl  diml  #  no  of  rows:  phase  encoding  direction 

header  -i  d2  dim2 
header  -i  tau  ti 

header  -i  naq  na  #  no  of  acquisitions 

#  header  -i  time_interval  date  #  debug  read  in  the  time  when  the  data 

#  was  acquired 

#  echo  "  time_interval  =  $time_interval"  #  debug 


# 

echo  "dO  = 

$d0" 

# 

debug 

# 

echo  "dl  = 

$dl" 

# 

debug 

# 

echo  "d2  = 

$d2  " 

# 

debug 

# 

get  -a 

# 

debug 

# 

echo  "tau 

=  $tau" 

# 

debug 

# 

echo  "naq 

=  $naq" 

# 

debug 

rm  -f  tl_$fn 
rm  -f  rawimdata 
rm  -f  rawimheader 
rm  -f  $fn.errmap 


#  view  1  1  $row  #  view  the  rowth  data  set  along  acquisition 

#  direction,  one  row  at  a  time 

#  printdata  -r  $col  1  >  %tmp  #  print  data  (real  and  imag)  into  file 

#  %tmp 

#  printdata  -r  $col  1  #  debug 

#  shellset  amp  =  %tmp  #  assign  the  data  point  (real  and  imag) 

#  to  ID  array  amp 

#  the  index  of  the  array  begins  from  1  and 

#  amp[0]  contains  a  space  character 

#  amp[l]  contains  the  real  part  and  amp [2] 

#  contains  the  imaginary  part  of  the  data 


#  Normalize  the  signal  intensity  for  multiple  acquisitions  to 

#  keep  correct  magnitude 

#  cal  -o  int  $amp[l]  $naq  div  #  int  =  amp[l]/naq 

#  echo  "  amp[l]  =  $amp[l]  int  =  $int  "  #  debug 

#  cal  -o  int  $int  100  mul 


#  Obtain  the  amplitude  at  the  logest  delay  as  the  initial  estimate  of  A 
header  -i  naq  na 

view  1  1  $row  #  view  the  rowth  block  FFT  transformed  data  set 

#  along  acquisition  direction 
printdata  -r  $col  1  >t.emp 
shellset  A_initial  =  temp 
cal  -o  A_initial  $A_initial [ 1]  $naq  div 
rm  -f  temp 

cal  -o  A_initiall  $A  initial [1]  Ine 


cal  -o  m  $n  1  sub  #  m  =  n  -  1 

echo  "  n  =  $n  "  #  debug 

header  -i  tau  ti 

cal  -o  tau  $tau  1000000  div  #  tau  =  tau  /  1000000 
view  1  1  $row  #  view  the  data  with  the  longest  delay, 

#  which  is  on  the  top  of  the  stack 

#  printdata  -r  $col  1  >  amp_longest_tmp  #  at  the  moment.  Then  read 

#  the  data  into 

amp_l o n g e  s  t_tmp 

#  shellset  amp_longest  ==  amp_longest_tmp  #  amp_longest_tmp  is  a  file 

#  rm  -f  amp_longest_tmp 


4  f  * 


#  Write  the  text  and  date,  into  input  file:  clear  rate  Isq  in 
rm  clear_rate_lsq_in 

@  col_tmp  =  $col  +  1 

cat  $time_interval  >!  time_temp  #  read  the  time  interval  and 

#  store  in  array  time_inter_temp 
shellset  time_inter_tempi  =  time_temp  #  new  code 

#  echo  "  col_tmp  =  $col_tmp  "  #  debug 

echo  "Lsqfit  Input  Data"  »clear_rate_lsq_in 

echo  "File  Base:  $fn  Voxel:  (  X  =  $col_tmp,  Y  =  $row  )" 
»clear_rate_lsq_in 

echo  "Fit  to  Clearance  Moded :  ln(Q(x))  =  ln{Q0)-  x/T  "  »clear_rate_lsq_in 

echo  "ClearanceModel  2  ln(Q0)-x/T  "  »clear_rate_lsq_in 

echo  "QO  $A_initiall  1  1/QO  "  »clear_rate_lsq_in 

echo  "T  700  1  x/ (T*T)  "  »clear_rate_lsq_in 

echo  "Time (mins)  In ( Signal_Intensity)  "  »clear_rate  Isq  in 

echo  "  $time_inter_temp [ 1]  $A_initiall  "  »clear_rate_lsq_in 


@  exl  =  1 
while ($exl  <=  $m) 

#  EXCHANGE  POSITIONS  OF  FILE  $EX1  TO  ZERO  (TOP)  OF  STACK 
@  index  =  $exl  +  1 

stackmv  $exl  0 

#  get  -a  #  debug 

#  READ  PIXEL  MAGNITUDE  INTO  SHELL  VARIABLE 

#  get  -a  #  debug 

view  J.  1  $row 

printdata  -r  $col  1  >  %tmp 
shellset  ampl  =  %tmp 
rm  -f  %tmp 

#  READ  HEADER  VARIABI.ES  FOR  T  FIT 

#  header  -i  tau  ti  #  tau  =  ti 

header  -i  naq  na  #  naq  =na 

cal  -o  ampl[l]  $ampl[l]  $naq  div  #  ampl[l]  =  ampl[l]/naq 

cal  -o  ampl  $ampl[l]  Ine  #  ampl  =  In (ampl [1] 

#  cal  -o  tau  $tau  1000000  div  #  tau  =  tau/1000000 

#  echo  "tau  =  $tau  Amplitude  =  $ampl[l]"  #  debug 

#  @  time_scale  =  $time:_interval  +  $time_scale 

echo  "  $time_inter_temp [ $index]  $ampl  "  »clear_rate_lsq_in 
stackmv  0  $exl 
@  exl++ 
end 

#  ENDWHILE 

cat  clear_rate_lsq_ir,L 
echo  "" 

Isqtool  clear_rate_l£ q_in 

echo  -n  "  Do  you  vrant  to  fit  another  voxel  (y/n)  ?  ->  " 
set  fit_again  =  $< 
echo  "" 
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•  t 


if  (  $fit_again  !=  y  &&  $fit_again  !=  Y  )  then 
echo  "" 

echo  "  Bye  Bye  !  " 
echo  "" 
exit ( 4 ) 
else 

@  time_scale  =  0 
@  fit_another_voxel  =1 
endif 


end 

#  end  of  first  while 

unset  n  exl  ex  fit_another_ voxel  col_tmp  time_interval 


50 


REPLY  TO 
ATTENTION  OF. 


DEPARTMENT  OF  THE  ARMY 

us  ARMY  MEDICAL  RESEARCH  AND  MATERIEL  COMMANC 
5C4  SCOTT  STREET 

FORT  DETRICK,  MARYLAND  21  <0?  Sil'  ' 


MCMR-RMI-S  (70-ly) 


26  Nov  02 


MEMORANDUM  FOR  Administrator,  Defense  Technical  Information 
Center  (DTIC-OCA) ,  8725  John  J.  Kingman  Road,  Fort  Belvoir, 
VA  22060-6218 

SUBJECT:  Request  Change  in  Distribution  Statement 


1.  The  U.S.  Army  Medical  Research  and  Materiel  Command  has 
reexamined  the  need  for  the  limitation  assigned  to  technical 
reports  written  for  this  Command.  Request  the  limited 
distribution  statement  for  the  enclosed  accession  numbers  be 
changed  to  "Approved  for  public  release;  distribution  unlimited. 
These  reports  should  be  released  to  the  National  Technical 
Information  Service. 


2.  Point  of  contact  for  this  request  is  Ms.  Kristin  Morrow  at 
DSN  343-7327  or  by  e-mail  at  Kristin.Morrow@deP.amedd.army.mil 


FOR  THE  COMMANDER: 


Enel 


% 

Chiidf  of  Staff  for 
ormatjJ.on  Management 


ADB263708 

ADB257291 

ADB262612 

ADB266082 

ADB282187 

ADB263424 

ADB267958 

ADB282194 

ADB261109 

ADB274630 

ADB244697 

ADB282244 

ADB265964 

ADB248605 

ADB278762 

ADB264450’ 

ADB279621' 

ADB261475 

ADB279568 

ADB262568- 

ADB266387 

ADB279633 

ADB266646 

ADB258871- 

ADB266038 

ADB258945 

ADB278624^ 


